Quantum damped oscillator II: 
Bateman's Hamiltonian vs. 2D Parabolic Potential 

Barrier. 

Dariusz Chrusciriski 

lO ■ Institute of Physics, Nicolaus Copernicus University 

ul. Grudziadzka 5/7, 87-100 Torufl, Poland 

CM 

£3 ■ Abstract 

We show that quantum Bateman's system which arises in the quantization of a damped 
harmonic oscillator is equivalent to a quantum problem with 2D parabolic potential bar- 
rier known also as 2D inverted isotropic oscillator. It turns out that this system displays 
the family of complex eigenvalues corresponding to the poles of analytical continuation of 
the resolvent operator to the complex energy plane. It is shown that this representation 
is more suitable than the hyperbolic one used recently by Blasone and Jizba. 
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In the previous paper [1] we have investigated a quantization of a ID damped harmonic 
oscillator defined by the following equation of motion 

X + 27X + KX = 0, (1-1) 



where 7 > denotes the damping constant. To quantize this system we follow an old 
observation of Bateman [2] and double the number of degrees of freedom, that is together 
with (1.1) we consider 

y - 2 iy + K y = , (1.2) 

bi ' 

i.e. an amplified oscillator. The detailed historical review of the Bateman idea may be found 
in [3]. For more recent papers see e.g. [4] and [5]. The enlarged system is a Hamiltonian one 
and it is governed by the following classical Bateman Hamiltonian: 

H(x,y,p x ,p y ) =p x p y -l(xp x -yp y ) + uj 2 xy , (1.3) 



where u = \J k — 7 2 } Now, performing a linear canonical transformation (x,y,p x ,p y ) — > 
(xi,x 2 ,pi,p 2 ): 

xi = —^= , pi = -y/ujy (1.4) 



X 2 = -\/wi , p 2 = --7= , (1-5) 



throughout the paper we shall consider the underdamped case, i.e. n > j 2 
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and applying a standard symmetric Weyl ordering one obtains the following quantum Hamil- 
tonian 

H = ujpAx — 7p0x, (1.6) 
where x = (xi,X2), P = (fii,P2) and we define two natural operations: 

1 2 

p A x = pxx.2 - P2X1 , p0x = x0p = - ^2(x k p k + p k x k ) . 

k=l 

Note, that [p A x, p x] =0. This operator was carefully analyzed in [1]. In particular it 
was shown that the family of complex eigenvalues 

h\0 = *4lf£> , (1-7) 

with 

E± = hwl±ih 1 {\l\+2n+l) , (1.8) 

found already by Feshbach and Tikochinsky [6], corresponds to the poles of the resolvent 
operator fi(H,z) = (H — z) —1 . Therefore, the corresponding generalized eigenvectors |f~H) 
may be interpreted as resonant states of the Bateman system. It shows that dissipation of 
energy is directly related to the presence of resonances. 

In the present paper we continue to study this system but in a different representation. 
Let us observe that performing the linear canonical transformation (x, p) — ► (u, v): 

x=^, P = ^, (1.9) 

one obtains for the Hamiltonian 

H = uir Au + H iho , (1.10) 

where 

£iho = ^(v 2 - 7 2 u 2 ) , (1.11) 

represents a Hamiltonian of a 2D isotropic inverted harmonic oscillator (iho) or, equivalently, 
a 2D potential barrier — 7 2 u 2 . Now, lov A u generates an SO(2) rotation on (u±, U2)-plane. 
Therefore, in the rotating frame the problem is described by the following Schrodinger equa- 
tion 

ihlptf = Hiho^ri , (1-12) 

where the rotating frame wave function ip^ = exp(iuj v A ut/K) tp. 

A ID inverted (or reversed) oscillator was studied by several authors in various contexts 
[7, 8, 9, 10, 11, 12, 13]. Recently, this system was studied in the context of dissipation in 
quantum mechanics and a detailed analysis of its resonant states was performed in [14]. The 
present paper is mostly devoted to analysis of a 2D iho. We find its energy eigenvectors and 
show that they are singular when one continues energy into complex plane. The complex 
poles correspond to resonant states of the 2D potential barrier [15, 16]. 

Finally, we analyze the Bateman system in the hyperbolic representation used recently in 
[5] by Blasone and Jizba. It turns out that this representation in not appropriate to describe 
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resonant states and hence the family of generalized complex eigenvalues found in [5] is not 
directly related to the spectral properties of the Bateman Hamiltonian. We stress that it does 
not prove that these representation are physically inequivalent. Clearly they are. Different 
representation lead to different mathematical realization which is connected with different 
functional spaces and different boundary conditions. These may lead to different analytical 
properties and hence some representation may display resonant states while others not. 

Prom the mathematical point of view the natural language to analyze the spectral proper- 
ties of Bateman's system is the so called rigged Hilbert space approach to quantum mechanics 
[17, 18, 19, 20]. We show (cf. Section 4) that there are two dense subspaces &± 6 L 2 (R 2 ) 
such that restriction of the unitary group U(t) = e ~ lHt / h to $± does no longer define a group 
but gives rise to two semigroups: U~(t) = U(t)\$_ defined for t > and U+(t) = U(t)\$ + 
defined for t < 0. It means that the quantum damped oscillator corresponds to the following 
Gel'fand triplets: 

$± C L 2 (R 2 U ) C $' ± , (1.13) 
and hence it serves as a simple example of Arno Bohm theory of resonances [20]. 

2 2D inverted oscillator and complex eigenvalues 
2.1 2D harmonic oscillator 

Let us briefly recall the spectral properties of the 2D harmonic oscillator (see e.g. [21, 22]): 

h 2 O 2 

tfho = -yA 2 + ^p 2 , (2.1) 



where the 2D Laplacian reads 



A -- — + -— + —— (22) 
2 dp 2 p dp p 2 dtp 2 ' 



and (p, (p) are standard polar coordinates on (u±, U2)-plane. The corresponding eigenvalue 
problem 

H bo il% = e%tfS, (2-3) 

is solved by 
where 

iltp 

$,(¥>) = -=, / = 0,±1,±2,... , (2.5) 

V Z7T 



and the radial functions 

Rni(p) = c nl (y/n/hp)W exp(-n P 2 /2h) 1 F 1 (-n,\i\ + i,n P 2 /h) , (2.6) 

where the normalization constant reads as follows 



_ ^2g/g (n + m „_ n 1 (27) 

C "'- mi V ~\ ' n - 0,1,2,... . (2.7) 
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Finally, the corresponding eigenvalues are given by the following formula 

e% = m(\l\ + 2n + l) . (2.8) 

Note, that using well known relation between confluent hypergeometric function \F\ and 
generalized Laguerre polynomials [24, 26] 

L ^- vllTi)Z% iFi{ - nMU (2 - 9) 

one may rewrite R n i alternatively as follows 



R nl (p) = ^2njh^ (n+'|/|)! (v/ ^ ,9)l " eM-n P 2 /2K)LM(n P 2 /h) . (2.10) 

It is evident that the family iff° is orthonormal 

(V#lA)=<W<^ , (2-11) 

and complete 

OO CO -. 

E E WpMtfaV* = -6(p- (t)S(<P ~ *0 , (2-12) 

n=0«=-oo ^ 

where ( | ) denotes the standard scalar product in the Hilbert space 

H = L 2 (R + ,pdp) ® L 2 ([0, 2n),dip) . (2.13) 

2.2 Scaling and complex eigenvalues 

Let us note that H^o defined in (1.11) corresponds to the Hamiltonian of the harmonic oscil- 
lator with purely imaginary frequency f2 = ±ry. The connection with a harmonic oscillator 
may be established by the following scaling operator 

V x -expQ v0uj , (2.14) 

with A 6 1. Using commutation relation [ufc,t)j] = ihSki, this operator may be rewritten as 
follows 

VA = e- iA exp(-iAp|^ , (2.15) 

and therefore it defines a complex dilation, i.e. the action of V\ on a function ip = ip(p, <p) is 
given by 

V x ip(p,<p) =e- iA V(e- a p^) • (2-16) 

In particular one easily finds: 

Vx H iho V' 1 = \ (-ft 2 A 2 - e- 4 * A 7 V ) • (2.17) 
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Therefore, for e 4tX = — 1, i.e. A = ±7r/4, one has 

Ir - 2 



Now, let us introduce 
that is 

It is evident that 
where 



V ±n/4 H iho V-* /4 = ±* (~y A 2 + ^ p 2 ) . (2.18) 

^ = ^/4*, (2.19) 

u± (p, = v/±? ^(V±ip, ¥>) • (2-20) 

^ihcu± =e±u± , (2.21) 



4 = ±^n? = ±^7(KI+2n + l) . (2.22) 

We stress that V\ is not unitary (for A € E) and hence in general Vj^ does not belong to H 
even for ip G H. In particular the generalized eigenvectors do not belong to H (the radial 
part R n i(\f±ip) is not an element from L 2 (E + , pdp)). 

Proposition 1 Two families of generalized eigenvectors satisfy the following properties: 
1. they are bi-orthonormal 

U nl(P' ¥>) U n'/'(^' ¥>) pdpcfy? = <W<5jj/ , (2.23) 

5. i/tey are bi- complete 

oo oo 1 

£ £ uS^) (p', ^') = -6{p- /)8(<P - tf/) . (2.24) 

n=0l=— oo 

The proof follows immediately from orthonormality and completness of oscillator eigenfunc- 
tions ?/>k°. 

3 Spectral properties of the Bateman Hamiltonian 

Now, we solve the corresponding spectral problem for the Bateman Hamiltonian (1.10). Note 
that H is bounded neither from below nor from above and hence its spectrum a(H) = 
(-co, oo). The corresponding generalized eigenvectors satisfy 

Hil> e , t = E e ^ x , (3.1) 

where I £ Z and e 6 E. Assuming the following factorized form of ip £i i 

1> e , l (p,<p) = R e , l (p)* l (tp) , (3.2) 

one has 

E £ + = uhl + e , (3.3) 
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with 

H iho R £ ,i = eRe,i . (3.4) 
The above equation rewritten in terms of (p, y?)-variables takes the following form 

and its solution reads as follows 

RsM = N e,l (\At7M 1 " exp(-i 1P 2 /2h) 1 F 1 (a, \l\ + l,i 1P 2 /K) , (3.6) 

with 

-i^l + l-jlj). (3.7) 
The normalization factor N £ j is chosen such that 

roc 

/ R~£p)R £ ,, l {p)pdp = 8{e-e') . (3.8) 
J o 

It turns out (see Appendix A) that 

Ne ' l = V^T Hrr(a) ' (3 " 9) 

with a defined in (3.7). 

Proposition 2 The family of generalized eigenvectors ip £j i satisfy the following properties: 

1. orthonormality 

ipe,i(p,v)ipe',i'(p,<p) pdpdtp = 5(e - £')8 lv , (3.10) 

2. completeness 

£ / de ^(p, ip) ^(p', <p') = -5(p- p')6(<p - if/) . (3.11) 

l=—oo 

Let us define another family of generalized energy eigenvectors 

Xe,i=?A,i , (3.12) 
where the anti-unitary operator T is defined as follows 



T^ l (p, ( p)=R Eil (p)^ l ( ( p) . (3.13) 

It easy to show that Bateman's Hamiltonian H is T-invariant 

THT* = H . (3.14) 

Moreover, if ip(t) = U(i)ipo, then Tip(t) = U(—t)(Tipo), which shows that T is a time 
reversal operator. Finally, Proposition 2 gives rise to the following spectral representation of 
the Bateman Hamiltonian 

-fOO poo +OQ /.qo 

H= / deE £ ^ l ){^ l \= Y / dsE £:l \ X£ ,i)(Xe,i\ , (3.15) 

l=—oo 

with E £j i defined in (3.3). 
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4 Analyticity, resolvent and resonances 



Now, we continue energy eigenfunctions tp £i i and Xe,i into complex e-plane. Note, that e- 
dependence enters R £ i via the normalization factor N £ j and the function ±Fi(a, \ l\ + 1, i^p 2 /K) 
(a is e-dependent, see (3.7)). It is well known (see e.g. [26]) that confluent hypergeometric 
function iF\ (a, b, z) defines a convergent series for all values of complex parameters a, b and 
z provided a / —n and b / — m, with m and n positive integers. Moreover, if a = —n and 
b ^ —m, then 1^1(0,6,2;) is a polynomial of degree n in z. In our case b = \l\ + 1 which 
is never negative and hence iF\(a, \l\ + l,i~/p 2 /h) is analytic in e. However, it is no longer 
true for the normalization constant N £ j given by (3.9). The T-function has simple poles at 
a = — n, with n = 0, 1, 2, . . ., which correspond to 

e = e n i = i 1 h{\l\+2n+l) , (4.1) 



on the complex e-plane. On the other hand the time-reversed function R £ i has simple poles 

at £ = £ n l = £nl- 

It is, therefore, natural to introduce two classes of functions that respect these analytical 
properties of ifj £t i and Xe,i- Recall [28] that a smooth function / = /(e) is in the Hardy class 
from above 7i\ (from below TL 2 _) if /(e) is a boundary value of an analytic function in the 
upper, i.e. Ime > (lower, i.e. Ime < 0) half complex e-plane vanishing faster than any 
power of e at the upper (lower) semi-circle |e| — ► 00. Define 



:= {0 G S(R 2 U ) I /(e) := (x £ ,i\cf>) ^^l}, (4.2) 



and 



$ + := [<P G S(R 2 U ) I /(e) := {^ e>l \<f>) £ 7i\ j , (4.3) 

where iS(M„) denotes the Schwartz space [29], i.e. the space of C°°(lRu) functions / = 
f(u\,U2) vanishing at infinity (|u| — > 00) faster than any polynomial. 
It is evident from (3.13) that 

$ + = T($_) . (4.4) 
The main result of this section consists in the following 
Theorem 1 For any function cfy^ G <£± one has 



^ + = E E U nl( U n^ + )^ (4-5) 

and 

f = EE Ki«i\<t>~)- (4-6) 

n=0 l=— 00 

For the proof see Appendix B. The above theorem implies the following spectral resolutions 
of the Hamiltonian: 



n=0 l=— 00 



00 00 



00 00 

H_ = H 

n=0 l=— 00 



E E E nlKl)Kl\ , (4-7) 
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and 



oo oo 



H + = H 



= E E E nlKl)(K 



n=0 l=—oo 



In the above formulae is given by (1.8). 

The same techniques may be applied for the resolvent operator 

1 



R(z,H) 



H-z 



(4.8) 



(4.9) 



One obtains 



oo „ 

R+(z,H) = e y 



de 



l=— oo 
oo oo 



E F i - z 



oo ^e,l 



\lp e ,l)(lpe,l 



= E E z Kl)( U nl\ > 



K" 1 — z 

n=0Z=-oo 



(4.10) 



on <I> + , and 



R-(z,H) 



de 



Eel - Z 



|Xe,l)<X. 



oo _ c 

S / 

( = — OO 
OO OO j 

E E E - _ \ U nl)( U nl\ > 
n=0/=-oo 



(4.11) 



on Hence, R+(z, H) has poles at z = E 1 ^, and R-(z, H) has poles at z = E 1 ^. As usual 
eigenvectors and u~ z corresponding to poles of the resolvent are interpreted as resonant 
states. Note, that the Cauchy integral formula implies 



^ := Kl)«l\ = ±f R+{z,H)dz, 

nl 



(4.12) 



where r+, is a clockwise closed curve that encircles the singularity z = E\. Similarly, 



(4.13) 



where T nl is an anti-clockwise closed curve that encircles the singularity z = E nl . One easily 
shows that 

Pnf P n>l> = 6nn>6iyP±, (4.14) 

and hence the spectral decompositions of (4.7) and (4.7) may be written as follows: 



oo oo 



^ = E E E ^i- 



(4.15) 



n=0 l=—oo 
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Finally, let us note, that restriction of the unitary group U(t) = e tHt / h to <&± no longer 
defines a group. It gives rise to two semigroups: 

(j_(t) := e-M-t/h . $ _ — „ $ _ ^ for t > , (4.16) 

and 

U + (t) ■= e ~ iA+t/h : — > $+ , for i < . (4.17) 
Using (4.15) and the formula for E n i one finds: 

OO OO 

<f,-(t) = U-(t)<p- = e-^ lt J2 e ~ l{2n+lll+1)tp ni > ( 4 - 18 ) 

£=— oo n=0 

for i > 0, and 

OO OO 

0+(t) = t>+(i)0 + = e"^* ^ e T(2n+|/|+l)tp+ ; (4>lg) 

i=— oo n=0 

for i < 0. We stress that <^ (4>f ) does belong to £ 2 (K 2 ) also for t < (t > 0). However, 
07 g (0+ g $+) only for t > (t < 0). This way the irreversibility enters the dynamics 
of the reversed oscillator by restricting it to the dense subspace $± of -L 2 (R 2 ). 

From the mathematical point of view the above construction gives rise to so called rigged 
Hilbert spaces (or Gel'fand triplets) [17, 18, 19, 20]: 

LcHcC, (4.20) 

and 

$ + CHCf' + , (4.21) 

where <&± denote dual spaces, i.e. linear functionals on <3?-|-. Note, that generalized eigen- 
vectors are not elements from TC. However, they do belong to The first triplet 
($>_,H, $'_) is corresponds to the evolution for t > 0, whereas the second one (& + ,H,& + ) 
corresponds to the evolution for t < 0. 



5 Bateman's system in hyperbolic representation 

In a recent paper [5] Blasone and Jizba used another representation. They transform Bate- 
man's Hamiltonian (1.3) into the following form 

H{yi,y 2 , w 1 ,w 2 ) = - (wf - wl) - j(yiw 2 + y 2 wi) + ^ 2 (yj - y\) , (5.1) 
with the new positions 

x+y x-y 



and new canonical momenta 



Px+Py Px-Py ,-„x 

. 1 = - 7r , ^ = ^r- (5 - 3) 
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Now, introducing hyperbolic coordinates (g,u): 



Hi = Q cosh-u , 1)2 = Q sinhn , (5-4) 

the canonical quantization leads to the following Hamiltonian defined on the Hilbert space 
H = L 2 (R + , gdg) <g> L 2 (R, du): 

H = H + H iho , (5.5) 



with 



H = -^D 2 + ^g 2 , (5.6) 



and the iho part 

H iho = i 7 h-^. (5.7) 
In the above formulae U2 denotes the 2D wave operator, that is 

□ - J^1 + IA_ ( 58 ) 

dy\ dy 2 dg 2 g dg g 2 du 2 

Clearly, in the (g, u) variables the formula for H-^ considerably simplifies and H^ represents 
the generator of SO(l,l) hyperbolic rotation on the (yi, ?/2)-plane. In this particular repre- 
sentation H\\ 10 defines a self-adjoint operator on L 2 (R, du). The corresponding eigen-problem 
is immediately solved 

Hfro$v = lftv®v , (5.9) 

with <& v (u) = e~ %uu /\/2n , and hence it reproduces the continuous spectrum of 2D iho 
<j(.Hiho) = (—00,00). However, there is a crucial difference between elliptic (p,ip) and (g,u) 
representations. The generalized eigenvectors & u may be analytically continued on the entire 
complex i/-plane. Therefore, the hyperbolic representation does not display the family of 
resonances corresponding to complex eigenvalues e n i defined in (4.1). Of course one may by 
hand fix the values of v to v = i(2n + \l\ + 1) but then the corresponding discrete family 
is neither bi-orthogonal nor bi-complete (cf. Proposition 2). 

To show how the complex eigenvalues of Blasone and Jizba [5] appear let us consider i^o 
defined in (5.6). Note, that Hq resembles 2D harmonic oscillator given by (2.1). There is, 
however, crucial difference between Hq and i^ho- The hyperbolic operator '—□2', contrary to 
the elliptic one '— A2', is not positively defined and hence it allows for negative eigenvalues. It 
is clear, since in the elliptic (p, ip) -representation Hq = iojhd v defines a self-adjoint operator 
on L 2 ([0, 27r), dtp) with purely discrete spectrum ujfd (I € Z). Now, the spectral analysis of 
the Bateman Hamiltonian represented by (5.5) is straightforward: 

Hip tv = £ eu ^ eu , (5.10) 

with 

£ eu = e + jhu , (5-11) 
and the following factorized form of ip €U : 

ip ev {g,u)=n ev {g)^ v {u) . (5.12) 
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The radial function TZ eu solves 



H K eu = eK eu , (5.13) 

and in analogy to (3.6) it is given by 

n ev (g) = N ev {^Jhg) iu exp(-LOQ 2 /2h) U(b, iv + 1, ug 2 /h) , (5.14) 

with 

b= 1 -(iu+l-4-) . (5.15) 



2 V hw, 

In (5.14) we have used instead of the standard confluent hypergeometric function \F\ so 
called Tricomi function U (see e.g. [27]). 2 It is defined by 

U(a, c, z) = r( ^~^ 1} i*i(a, c, z) + ^ ~ ^ z x ' c lFl (a-c+ 1,2- c, z) . (5.16) 

A Tricomi function U(a, c, z) is an analytical function of its arguments and for a = —n 
(n = 0, 1, 2, . . .) it defines a polynomial of order n in z: 

U(-n,a + l,z) = {-l) n n\L«{z) . (5.17) 

Moreover, using the following property of U (an analog of (B.5) for iFi) 

U(a, c, z) = z 1 - E7(l + a - c, 2 - c, z) , (5.18) 

one obtains 

/•oo fc POO 

/ K ev (Q)K el/ ( e ) gdg = — \N el/ \ 2 / z lv e' z U 2 (b, iv + 1, z) dz , (5.19) 
Jo ZuJ Jo 

with z = ujg 2 /h. Now for b = —n, lZ ei/ belongs to the Hilbert space L 2 (M + , gdg). It implies 

e = hw(2n+ 1 + iv) . (5.20) 

and hence it reproduces discrete spectrum l fkv x integer' iff iv = I = 0,±1,±2, .... Now, 
using (5.17), (5.19) and 

f°° e- z z a LZ(z)L^(z) dz = - r(n + a + 1) 5 nm , (5.21) 
Jo n - 

with a > — 1, one obtains the following family 1Z n i G L 2 (R + , gdg): 3 

^ = ^/ n!r (n 1) (VWM* exp(-^^ 2 /2/ l ) L l n (ug 2 /h) , (5.22) 

2 Actually, in [27] (and also in [23]) this function is denoted by G. We follow the notation of Abramowitz 
and Stegun [26]. 

3 There is a difference in normalization factor in formulae (37) in [5]. It follows from slightly different 
definition of L" . 



K 



nl \ 
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with n = 0, 1, 2, . . . , and I = 0, 1, 2, . . . , satisfying 



/ TZ n i(g)TZ n n(g) gdg = 6 nn > . (5.23) 
Jo 

We stress that the family lZ n i is defined for I > only (otherwise it can not be normalized!). 
Finally, defining 

cf> nl (g,u) = -l=H nl {g)e~ ul , (5.24) 

V27T 



one has 

H4>nl = £nl 4>nl , (5-25) 

with 

£ nl = frw(2n + I + 1) - itvy I . (5.26) 

There is, however, crucial difference between families u^(/?, (p) and 4> n i(g, u). The family |u^) 
corresponds to the poles of ij) £i i from (3.2). No such correspondence holds for \<j> n i) and 
from (5.12). In particular there is no analog of Theorem 1 for \4> n i). Moreover, £ n i contrary to 
E n i from (1.8) does not fit the formula for complex eigenvalues of Feshbach and Tikochinsky 
[6] (see detailed discussion in [1]). It defines simply another family which is however not 
directly related to the spectral properties of the Bateman Hamiltonian. 

Appendix A 



To compute N £ j in (3.6) let us analyze the quantity I £ = Jq°° R £t i(p) R £i i(p) pdp. Clearly, this 
integral diverges (I £ = <5(0)), however, its structure enables one the calculation of N £ j. One 
has 

1 f'°° 

h = ^\N £ ,i\ 2 z\ l \ l F l {a,\l\ + l,iz) l F l {a,\l\ + l,-iz)dz , (A.l) 
2 7 Jo 

where we defined z = jp 2 . Now, the integral in (A.l) belongs to the general class 

poo 

J= e~ Xz z^ 1 1 F 1 (a,p,kz) 1 F 1 (a',p,k'z)dz , (A.2) 
Jo 

given by the following formula (see Appendix f in [23]): 

J = r( M ) A Q - Q '^(A - ky a (X - k')- a ' 2F1 (a, a', p- . (A.3) 
Using the above formula with A = 0, p = \l\ + 1, a = a, a' = a and k = —k! = i one finds 

h = ^ |Ay 2 (-r a Hp2fi(a,a, \l\ + 1; 1) . (A.4) 

Finally, noting that 

P , fl(V . n r( 7 )r( 7 -«-/?) 

2 F 1 (a, (J, 7; 1) = — — , 

r(7 - a)r(7 - 0) 
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one has 



4 = || Ne ,, PHr (_,)-. _M_ 



Therefore, comparing (A. 5) with I £ = 5(0) one finds 



7 



7T 



H) a r(a), 



(A.5) 



(A.6) 



which proves (3.9). 

Appendix B 

Due to the Gel'fand-Maurin spectral theorem [17, 18] an arbitrary function <p + £ <£> + may be 
decomposed with respect to the basis ip £ i 



°° />0O 



(B.l) 



/=— oo 



Now, since ( ip £) i\4> + ) £ we may close the integration contour along the upper semi-circle 
|e| — > oo. Applying the Residue Theorem one obtains 



-oo ra=0 



(B.2) 



Using the well known formula for the residuum 



ResT(a) 



(-l) r 



a=-n n\ 



one obtains 



Res V', 



6,1 



[n 



+ 1/1)! 



e=e nl ^i2n+\l\+l \ 2ttH ]j n\\l\\ nl 

Moreover, the analytical function ip £i i computed at e = e n i reads: 



(B.3) 



(B.4) 



'el 



.n+|i|+l _1 



^E(y/-i>y/hp)W exp(z 7/ 9 2 /2^)i^i(n + K| + l,|Z| + l,-HP 2 A)^M • 



Due to the well known relation [24, 25, 26] 

iF 1 (a, b, z) = e 2 iFi(6- a,b,-z) , 
one finds 



and hence the formula (4.5) follows. In a similar way one shows (4.6). 



(B.5) 
(B.6) 
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